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A special class of transmission lines is considered, in which the modes 
decompose into two noninter acting sets. Both a single transmission line 
with constant characteristic impedance and variable propagation factor, 
and two transmission lines with equal propagation factors and variable 
coupling, in which the forward modes do not interact with the backward 
modes, are investigated. Exact expressions are obtained for the reflection 
and transmission coefficients when a section of such a transmission system 
connects two semi-infinite transmission systems consisting of constant 
impedance and admittance lines. These results hold for arbitrarily varying 
propagation factors and coupling; and while they are of independent 
interest in the case of deterministic variations, we make an application of 
them here in the case of stochastic variations. 

Exact results are obtained for the ensemble averages of the transmission 
coefficient and transmitted power, and their variances, for the inserted 
section of single line, when the variable propagation factor is a random 
function involving either a Gaussian process or the random telegraph 
process. Asymptotic results are also obtained in the general case of weak 
fluctuations and long inserted sections. Analogous results may be obtained 
for the inserted section of two lines when they are randomly coupled, and 
the results are given in the case of matched lines, for which no reflections 
occur. Finally, some of the time domain statistics for lossless lines are 
considered, and expressions are derived for the ensemble averages of the 
transmitted pulse, due to pulses incident on the inserted section. 

I. INTRODUCTION 

This paper deals with a special class of the generalized equations of 
telegraphy. The starting point is the following simple observation. 
Consider the telegraphist equations in the frequency domain for a 
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single transmission line 

| = -<y(x)v(x), (l) 

where V(x) and d(x) are the time Fourier transforms of the voltage and 
current in the line, and d(x) and •yfx) are the impedance and admittance, 
respectively. Then, if the characteristic impedance 

k = Vmfim (2) 

is a constant independent of x, it is simple to show that (1) has the two 
fundamental solutions 

V + (x) = exp I- J W dt\ , 6 + (x) = jj-V + (x), 

V-(x) = exp || o " r© #} , 6-{x) - -j'O-Cr), (3) 

where 



T(x) - V^(i). (4) 

This simple result has several interesting consequences. Note that 
1) + (x) and 0+(x) describe a wave moving to the right and V-(x) and 
B-(x) .describe a wave moving to the left (the time factor is assumed 
to be «'"'). Furthermore, the wave moving to the right does not induce 
a reflected wave moving to the left, and vice versa, except possibly at 
the beginning or termination of the line. 

This decomposition into noninteracting right and left moving waves 
suggests that a similar decomposition may exist in the case of n coupled 
transmission lines. This is the subject of Section II. It is shown there, 
that when a condition analogous to (2) is satisfied by the impedance 
and admittance matrices of the system, the system of 2n equations can 
be decomposed into two noninteracting sets of n equations. Under some 
circumstances, the n fundamental solutions of one system correspond 
to waves moving to the right; and the n fundamental solutions of the 
other system correspond to the left moving waves. In particular, we show 
that the model of two interacting waves recently studied by Rowe and 
Young 1,2 corresponds to just such a decomposition. 

In addition, the form of the solution (3) for the single transmission 
line makes these models particularly convenient to study when d(x) 
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and 'y(x) are stochastic processes. In Section III, we consider a finite 
section of length / of lossless transmission line with variable impedance 
and admittance satisfying condition (2). This section of line connects 
two semi-infinite lossless transmission lines having constant inductances 
and capacitances. We derive various expressions for the reflection and 
transmission coefficients of the inserted section. Using these expressions, 
in Section IV, we calculate exact expressions for the ensemble average of 
T(o), I), the transmission coefficient, and of | T(a), I) | 2 and | T(u, I) | 4 , 
when d(x) and ^(x) are particular stochastic processes. Asymptotic 
results are obtained in the general case of weak fluctuations in d (x) and 
'y(x), and long inserted sections. These results are based on a limit 
theorem of Khas'minskii, 3 and the details are given in the Appendix. 

In Section V, we study a more complicated model consisting of two 
coupled transmission lines. Here the model involves a finite length, /, of 
two transmission lines having equal propagation constants, with variable 
inductive and capacitive coupling, connecting two semi-infinite, constant 
impedance and admittance lines. The self-impedance and admittance 
of the inserted section are also constant, and the semi-infinite lines are 
uncoupled. We derive expressions for the reflection and transmission 
matrices of the inserted section. Exact results may be obtained for the 
ensemble averages of the elements of the reflection and transmission 
matrices in the case of random coupling, for particular stochastic 
processes, using the results of Section IV. Since the results are quite 
lengthy in the general case, we give them only in the case of matched 
lines, so that no reflections occur. 

Finally, in Section VI, we consider some of the time domain statistics 
of our models. Exact expressions are derived for the ensemble average 
of the transmitted pulse, due to a pulse incident on the inserted section 
of single transmission line with random inductance and capacitance. It 
is of interest to note that if the fluctuations in the propagation factor 
are described by a Gaussian process, the transmitted wave violates 
causality. This is not the case when the fluctuations are described by 
the random telegraph process. Analogous results are obtained for the 
transmitted pulses, due to pulses incident on the inserted section of two 
randomly coupled transmission lines, in the case of matched, lossless 
lines. 



II. CLASS OF TRANSMISSION LINES 

The generalized equations of telegraphy are the starting point of 
this paper: 
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j- V(a>, t) = -R(x)I(x, t) - L(x) £ t I(fi, t), (5a) 

jj I(x, = - G(x)V(x, t) - C(x) f t V(x, 0- (5b) 

These equations are typically used to describe the time and space 
variations of the current and voltage along n coupled transmission 
lines. 4 In this case, V and I are column vectors whose elements V p (x, t) 
and I v (x, t), p = 1, 2, • • • , n, are, respectively, the voltage of the pth line 
relative to some fixed voltage, and the current in the pth line. R, L, G, 
and C are n X n matrices, the resistance, inductance, conductance, and 
capacitance respectively, which typically are functions of the distance 
x along the lines. 

We will, for the most part, find it convenient to work in the frequency 
plane, and so we introduce the Fourier transforms 



V(x, o>) = -^=- 1" V(x, t)e- iat dco, 



(6a) 
(6b) 



Then V and I satisfy the equations 



where 

Z{x) = R(x) + joL(z), Y(a:) = G(x) + juC(x) (8) 

are the impedance and admittance matrices, respectively. It should be 
noted that the frequency domain equations appear in other contexts, 8 
but there, the frequency dependance of Z(x) and Y(x) is generally more 
complicated. The remainder of this section is devoted to some general 
properties of the frequency domain eqs. (7) . 
It follows from (7) that 

4~ (VI* + V*'l) = -I'CZ'Cb) + Z*(z)]l* - V'[Y'(.r) + Y*(z)]V*, (9) 
ax 

where t denotes transpose and * denotes complex conjugate. Hence, 
for lossless lines 

Z'(x) + Z*Cr) = 0, Y'(x) + Y*(x) = 0. (10) 
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We wish to consider the class of coupled transmission lines for which 

Z(x) = KY(x)K, (11) 

where K is a nonsingular constant matrix. Note that, in the lossless case, 
the first condition in (10) follows from the second if K' = K*. When (11) 
is satisfied, the solutions of (7) may be split into two groups, namely 



and 



v = K, > i - - KY ^ V > 



v - - Kl ' £ - KY ^ V - 



(12) 



(13) 



If the lines are appropriately matched at both ends, then either one set 
of solutions or the other occurs, and reflections are avoided. Since K is 
constant, this matching is independent of the length of the lines. 
As a particular example, let n = 2 and 



K = 



K x 
.0 K 2 



Y = 



I\ -jc(x) 

-jc(x) 



.0W K 2 _ 



(14) 



Thus, from (11), 

z = [" M, -m&Mx) 

-JiK.KMx) T 2 K 2 

Corresponding to (12), we have 



(15) 



dx 

dV 2 
dx 



+ r x Vt = Mx)(jr]v 2 , 
+ T 2 V 2 = j&j^Jvx ■ 



(16) 



The substitutions D, = K\l 1 , V 2 = K\l 2 lead to the equations for two 
coupled modes traveling in the same direction, which have been con- 
sidered previously. 12 We remark that if we choose K = [£« ° K J instead, 
then we are led to equations for two modes traveling in opposite direc- 
tions. 

Next, we consider a particular class of transmission lines, satisfying 
(11), for which 



KY(*) = r(»)I - jc(x)A, 



(17) 
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where I is the unit matrix of order ??, and A is a constant matrix. Then 
there are solutions of (12) and (13) of the form 

V = b exp |=F j" [r(.r) - j\c(x)] dxj = ±KI, (18) 

where b is a constant vector satisfying 

(XI - A)b = 0. (19) 

The eigenvalues X are given by | XI — A | = 0. 

We will assume that e(z)A describes only the coupling between lines, 
so that A has diagonal elements equal to zero. The case n = 1 (for which 
A = 0, and X = is the only eigenvalue) corresponds to the well-known 
case of a single line with constant characteristic impedance and variable 
propagation factor. This case is considered further in subsequent 
sections. The case n = 2 corresponds to n transmission lines with 
identical propagation factors and variable coupling. Such a situation 
might arise in the consideration of n twisted pairs in a cable, although 
the relationship (11) is not too realistic. In Section V, we consider the 
case n = 2 corresponding to I\ = Y 2 in (14), so that (17) holds. 

III. SINGLE TRANSMISSION LINE 

In this and the following section, we study in some detail the following 
example. Consider an infinitely long, lossless, single transmission line 
in = 1 in the classification of Section II) which for x < has the 
constant impedance and admittance d = ji»L , ^o = jtaC , for x > I 
has the constant impedance and admittance 3, = jcoLj and ^ = juCi , 
while the central section < x < I has the variable impedance and 
admittance d(x) = jwL(x), <y(x) = jo>C(x). A wave traveling to the right 
in the region x < will be partially reflected and partially transmitted 
on striking the central region in < x < I. We study the transmitted 
wave under the assumption that 2>(x) and 'y (x) satisfy condition (11), i.e., 

K = Vd(x)/^(x) = VL(x)/C(x) (20) 

is a positive constant independent of x. 

Although this case is probably hard to realize physically, it is never- 
theless of considerable interest, since it is mathematically simple enough 
so that many interesting questions about it can be answered. 

Let 

r Q = Vd a *y a = juVL a C a = jay a , 

K a = VdJ%, = VLJC~ a , « = 0,1. (21) 
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Then in x ^ and x ^ I, we can solve eqs. (7) simply. In x ^ 0, we 
have a solution 

V(x) = e~ roX + R(a, l)e T ° z , 

${v) = IT ( e ~ r ° z ~ R(«>> l ) eT " x )- (22) 

■/v 

This represents a plane wave e~ v ° z moving to the right and a reflected 
wave R(co, l)e r ° z moving to the left, where R(u, I), the reflection coeffi- 
cient, is a function of co and /. Similarly, in x ^ I, the solution can be 
written 

V(x) = 7>, l)e- r,(z - n = K t <l(x), (23) 

representing a transmitted plane wave moving to the right, where 
T(u, I) is the transmission coefficient. 
In ;£ x ^ I we define the propagation factor 

T(x) = VWtfx) - jo>KC(x) = jco U&- (24) 

Then, from (18), we can write the general solution of (7) in ^ x ^ I as 
V(x) = AV + (x) + BV.(x), 

6{x) = ± (AV + (x) - BV.(x)), (25) 

when A and B are constants and 

V ± (x) = exp|=F£r(0^}- (26) 

We now have a solution depending on four unknown constants 
A, B, R and T which can be determined from the condition that V(x) 
and d(x) must be continuous at x = and x = I. The resulting four 
linear equations are easily solved and yield for the reflection and 
transmission coefficients 

m n _ (K - K )(K + K,)VM) - (K + K a )(K - K,)V + (l) m , 
«l«, V {K + Kn)(K + Kl)v _ { i) - ( K - K )(K - K t )V + (l) ' K "> 

ATJJf 

T( - U3 ' l) = {K + K )(K + K,)V_(l) - (K - K )(K - K,)V + (l)' (28) 

Notice that if K - K = K t , then R(w, I) = 0. 

We confine our study to the transmitted wave, although the reflected 
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wave can be studied equally well by the techniques we employ. In 
particular, there is the easily proved energy conservation relationship: 



i = |j?(«, or + f; \n<*, 01 



(29) 



Before proceeding, let us assume instead that a section of line 
^ x ^ I is driven by a voltage source in series with an impedance Z 
and that the line is terminated in an impedance Zi , as shown in Fig. 1. 
Then, V(x) and d{x) are still given by (25), but the boundary conditions, 
from which A and B are determined, are now 

-0(0) + Z o d(0) = E n , 

V(l) = Ztd®. (30) 

Then, it is easily shown that 

2KZ l E 



V(Q = 7i 



(3D 



(K + Z )(K + Z,yO-(J) -(K- Z )(K - Z t )VS) 

so the transfer impedance in this formulation is essentially identical 
with the transmission coefficient in the first formulation. We shall 
continue to use the first formulation. 

We now further specialize the model, and let 

L(x) = L(l + *N(x)), (32) 

where we assume that ^ e ^ 1 is a dimensionless constant, L > 
is a constant with the dimension of inductance, and N(x) is a (dimen- 
sionless) stochastic process with zero mean. It follows that L is the 
stochastic mean of L(x), 

L = (L{x)). (33) 

The symbol ( ) will be used throughout to denote the stochastic mean. 




"UiD.^cc) 



z* 



Fig. 1 — Diagram of transmission line circuit driven by source E at x =0 with 
internal impedance Z and terminated at x = I by impedance Zu 
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It follows from (20) that 

C(x) = C(l + eN(x)), K = ^- (34) 

Further define 

7 - VCL, (35) 

so that 

T(x) = juy(X + eN(x)). (36) 

In addition, we let 

m = f N® «i (37) 

■Jo 

and 

(g - g,)(g - K t ) = 4KK t 

9 (K + K )(K +K t )> (K + K )(K + Ki) ' w 

Then we can write expression (28) for T(ai, I) as 

r(«, i) = xr ,<,TC, **' cl>) [i - p8 -w-Tci + .«i»j-i i (39) 

We can now derive some series expansions of T(u), I) and some of its 
powers which will prove useful in the next section. Since K, K and K t 
are positive, 

I P I < 1, (40) 

hence we have the geometric series expansion 

r(«, = X E p r e -l*r+l)i*r<l + *Hm (41) 

r-0 

Next let 

f(<p) = X 2 /[(l - pT")(l ~ pe")l (42) 

Then /(<p) has the partial fraction expansion 

and hence /(<p) has the series expansion 

m = r^ E p ln, e'". (44) 

■l P n--oo 



(45) 
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Then, it follows that 

m*, or = i^m t |l .H - pri-»<.*t«» < 

i- P r = -oo 

We finally need to examine double sums of the form 

S= E E P ,r,+U, *(r + «)- (46) 

If we make the change of summation variables u = r + s, t = r, we can 
write S as 

S = E h(u)m, (47) 

U-— 00 

where 

m = i, P " i+ "- ui . (48) 

I--00 

The series $(u) can be evaluated simply, and so we obtain 

s = E (m + f ± ^V" , *W- (49) 

„ri L' i - p J 

If we now square expression (45) for | ^(aj, I) | 2 , we obtain a series of 
the form (46), and so it follows that 

in., or = (T^ t {w + l±^}v"^'"" + -«». (so) 

IV. FREQUENCY DOMAIN STATISTICS 

In this section we study some of the frequency domain statistics of 
the model described in Section III. In particular, we obtain expressions 
for the stochastic average of T(u, I) and | T(u, I) \ 2 and examine the 
standard deviation of these quantities. 

It is clear from (41), and (45) and (50) that the problem of calculating 
(T(u, 0), < | T(ta, I) | 2 ) and ( | T(u, I) | 4 ) has been reduced to the problem 
of calculating (e' aeil) ). 

Consider first the case where N(x) is a zero-mean, Gaussian random 
process. Then 6(1) is a zero-mean, Gaussian random variable with 
variance 6 



'(J) = f f (N(x)N(y)) dx dy, (51) 
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and 

<a*- #c "> = — 7= f «'«•-••'»•■ dfl = e -i*"\ (52 ) 

<r\/27r -'-co 

Consequently, 

Mi(», I) = (T( u , 0) - X E P' r 'e(2r + 1, a,, Z), (53) 



r-0 



Af a («, = <|r(«, 0|"> = =-^-| £ p lrl e(2r,co, Z), (54) 

1 — P r=-oo 

and 

M 4 (co, I) = (|T(«, 0|*> 

= (1 _ pV r S {m + r^v} plr ' e(2r ' w ' * } ' (55) 

where 

e(r, co, - exp { -jroryl - ±rV 7 2 eV(Z)}. (56) 

If N(x) is white noise, then (N(x)N(y)) = D 5(x — y) and 

cr 2 (0 = D Z , (57) 

where D is a constant having dimensions of length. If N(x) is a wide- 
sense stationary Gaussian process, 

(N(x)N(y)) = g(\x-y |), (58) 

with continuous gr($), then 

a 2 (l) = 2 [ (I- Qg(Q #. (59) 

In particular, if g(£) = e~ 2H , then 

• p »-ir , - i Tf =r ]- <w 

Since o-(O) = in all these cases, e(r, to, 0) = 1, and so 

Jf.(« f 0) = [jzrzj ' » - 1. 2, 4. (61) 

In many cases, such as (57) and (60), c 2 (Z) tends monotonely to » as 
Z — » oo. In these cases we obtain the asymptotics of the moments as 
Z — > oo simply. For o)yea(l) > 1, we have 
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M 1 (o>,l)ttU- iayl - iu ' y """ U) , (62) 

M t (», I) « f 3 7 {1 + 2 P cos (2a>yl)e- 2 "'->"'" a) } , (63) 

M 4 (co, « (1 ** p y {1 + P 2 + 4p cos (2o>yl)e- 2 »"'""° ll) } . (64) 

The mean-amplitude transmission coefficient M,(w, Z) decays exponen- 
tially to zero as I — > °° , while the mean-power transmission coefficient 
tends to the limit 

X 2 4KK 2 (65) 



1 - p 2 " (Xo + X,)(K" + K Kd 

This can be explained qualitatively by noting that, when the transmitted 
amplitudes are averaged over the ensemble, cancellation can take place, 
while the transmitted powers all have the same sign and so no cancella- 
tion can take place on averaging. 

It is easily seen from (62) through (64) that the ratio of the standard 
deviation to the mean of T(a>, Z) is 

Z,(«,PW /} , «*-*"*•">, (66) 

VI — p 

while the ratio of the standard deviation to the mean of | T(to, Z)| 2 is 

2 3 («, Z) « ^^ [1 " P cos (2o> 7 0e-'" ,r ' , "' a) ]- (67) 

For the examples (57) and (60), S x («, Z) -* °o exponentially as Z -> oo, 
while S 2 (w, Z) tends to the limit 

"27" I K - K | | K - g, | 



(68) 



\1 - P 2 " V2K(K + K,)(tf 2 + JCo^i) 
To get some feel for the numbers, we note that for 

!*§ -|^ 2,0 S^£ 0.168. 

However, if if and if , differ too much from K, this ratio becomes much 
larger than 1. 

As a second example, consider the case where N(x) is the random 
telegraph process. 7 It is an ensemble of square wave functions {N(x)}, 
such that each sample function N(x) can assume only the values ±1. 
For fixed x, a sample function chosen at random will equal +1 or — 1 
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with equal probability. The probability p(n, x) of a given sample 
function changing sign n times in an interval of length x is given by the 
Poisson process 



(hx)* ... 



p(n,x) = K -^-e- bx , (n = 0, 1,2, •••)• (69) 



This process has zero mean and correlation function 

(N(x)N(y)) = e - 26| -«". (70) 

The probability density function for the integral of the random tele- 
graph process has been derived by McFadden, 8 and in our notation 

P(e,i) = K 6 Tw- 0) + 6(1+ 0) 

+ b{i [b(? - a 2 )*] + u ffiSfjfi} B Q - W + *>] » ( 71 > 

where H(u) is the Heaviside function 

H(u) = 1, u > 0, 

ff(tt) = 0, u < 0. (72) 

It follows that 9 

(«'-•<») = e - 6, rco8h {(6 2 -a 2 )*/} 

+ (b»- g y Binh ^ ~ « 2 )' l ' (73) 

If we define 
/(r, to, = exp {— jruyl — bl\ 

Tn^oK \<tf J 2 2 Mn _l b8mh {(b 2 - rVyVyzn ».. 

• cosh {(6 — r w 7 e )'i) H „ 2 _ jn »h| J » (74) 

then the expressions for (T(u, I)), (| T( W| Z)| 8 ) and (| T(a>, Z)| 4 > for the 
random telegraph case can be obtained from (53) through (55) on 
on replacing e(r, w, l) by /(r, «, I). 

If e > is small enough so that for a given positive integer r 

v = (er u>Y/b) 2 « 1, rfbl « 1, 

then for ^ r ^ r , 

i(r, co, I) = exp { -jraryl - frWy 2 eH/b\[l + 0(u) + 0(n'6l)]. (75) 
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If furthermore bl » 1, then from (60) for the Gaussian with correlation 
function (70), we have 

Hence, for these two cases, for ±? r ^ r , 

^4 = [1 + 0(„) + Otfbl)]. (76) 

e(r, w, I) 

Note that both e(r, w, I) and /(r, w, Z) are exponentially small for r > r 
if v bl is moderately large. Also, from (74), as I -* w the first-order 
moment of T(co, Z) tends to zero, and the second moment and its standard 
deviation tend to the same limits as in the Gaussian case. 

We now consider a general case of weak fluctuations in the inductance 
and capacitance, and long sections of line, so that < e « 1 in (32), 
and I = A/e 2 . It is assumed that N(x) is a bounded, zero-mean, wide- 
sense stationary stochastic process, with correlation function given 
by (58). An application is made in the Appendix of a limit theorem due 
to Khas'minskii, 3 in order to determine the behavior of 



E 0(A/e 2 ) = e / N(x) dx, 

Jo 



(77) 



for A bounded and e — > 0. 

If the stochastic process N(x) satisfies a certain strong mixing condi- 
tion, 3 then it is found that the process ed(A/e 2 ) converges weakly to a 
Markov diffusion process ®(A), with probability density function 

p< e - A > - &s? exp [ifi ■ (78) 



where 



(79) 



a = 2 Urn [^ j (X - z)g(z) cfej , 

and g is given by (58). If, as we assume, 

lim [^ f*zg(z) dfe] = 0, (80) 

then 

a = 2 r g(z) dz. (81) 

J a 
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It follows from (78) that 

(exp [;r«760(A/e 2 )]> PH exp (-±drWy 2 A). (82) 

Hence, for < £ « 1 and I = A/e 2 , the asymptotic expressions for 
(T(a, I)), <| T(ju, l)\ 2 ) and (| T(w, Z)| 4 ) can be obtained from (53) through 
(55) on replacing e(r, w, I) by 

h(r, w, I) = exp [-jruyl - £rV 7 YdZ}, (83) 

where d is given by (81). Fc* the random telegraph case g(z) = 
exp ( — 262), so that d = 1/6, anu the consistency of (75) with the above 
result is noted. 

V. TWO COUPLED TRANSMISSION LINES 

We consider here the case of two coupled transmission lines described 
in eqs. (14) through (16), but with identical propagation constants 
so that 

r, = r = r 2 . (84) 

Then V is independent of x in (17), and 

. r (KjK 2 f 

A = 

_{K S /K$ . 

Thus, the eigenvalues of A are X = ±1. Let 



(85) 



&0 = f c(y)dy, (86) 

Jo 

it being assumed that c{x) is real, and that K X K 2 > 0. Then, from (14), 
(18), and (19), the general solution of (7), subject to (11), for this case 
may be written in the form 

1), = K\[e m '\Be- T * + Ce Tx ) + (T ii{x) {Ae" Ts + De Tx )], (87) 

V 2 = K\[e m '\Be- T ' - Ce Vz ) - »-" w ^- r " - De T ')] t (88) 

6l = K:\e> iM (Be- Tz - Ce Vl ) + e' iiu) {Ae~ v ' - De T ')], (89) 

3 2 m K?[e m '\Be- Tx + Ce Vl ) - «- ,lw (4T r- + De Tl )}. (90) 

We suppose that the coupled lines extend from x = to x = I, so 
that (87) through (90) hold for g x g I. For x < and x > I we 
suppose that the transmission system consists of uncoupled lines with 
constant propagation constants r p(l and T pt , and constant characteristic 
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impedances K pn and K pl , respectively, (p - 1, 2). Then, for an incoming 
wave on line 1, 

Vl = ( e - r »' + R/" 1 ), V 2 = R 2 e r " x , (91) 

9l = 1 L (6 -r». _.«/-), <r 2 = ^a p -, (92) 

for z ^ 0, and 

^ m r/ r lI ««-«) ) ^ = r ^- r " to - ,, l (93) 

^-l-La-™-", 4 2 = £ 2 -e- r '"*-", (94) 

for x ^ I. 

The boundary conditions are that 1) p and 0, , (p = 1, 2), must be 
continuous at z = and x = Z. The calculation of reflection and trans- 
mission coefficients is tedious, but straightforward, so we omit the 
details and merely state the result. Let 



K p0 _ K vi 



(p-1,2), (95) 

A„ /v p 

and 



. = (Kt/Ktf, x = K0- W 

Also, define 

A = [(1 + m.)(1 + M 2 )(l + rO(l + *2)e 2r ' ~ 2(/x, - ft)fa - fO 

+ (1 - Ml)(l ~ mOU - "0(1 - "2K 2r ' 

- 2(1 - m.M 2 )(1 - "i"0 cos 2x]. (97) 
Then it is found that 

AT, = ^[(1 + ih)(l + "0^' - (1 - m 2 )(1 - " 2 >" r '] coax, (98) 

AT 2 = 4jw 2 [(l + m 2 )(1 + *y + (1 - mOQ - ^>~ ri ] sinx, (99) 

as, = [(1 - m.)(i + m 2 )(i + ".)(i + *y vl + 20*1 + feXn - "«) 

+ (l + m0(1 - m0(i - "0(1 - " 2 >" 2ri 

- 2(1 + m.m 2 )(1 - *xv%) cos 2x1 (100) 



and 



AR* = -4i«/i 2 (l - "i"») sin 2 X . (101) 
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We remark that the reflection and transmission coefficients correspond- 
ing to an incoming wave on line 2 may be obtained by appropriate 
interchange of subscripts. 

Since the coupling function c(x) is real, x is also real, from (86) 
and (96). If k is an integer, then T a = and R 2 = for x = kir, and 
T t = and R 2 = for x = (k + h)ir. The oscillatory behavior of 
7\ and T 2 has been found earlier by Foschini, 10 in the case of matched 
lines with equal characteristic impedances, i.e., /*„ = 1 = v p , (p = 1, 2), 
and k = 1. In this case there are no reflections. 

In the general case, the expressions for the transmission and reflection 
coefficients 7\ , T 2 and R t , R 2 may be expanded in Fourier series in x- 
Thus, in the case of random coupling between the lines, the problem of 
calculating the expectations of the transmission and reflection coefficients 
reduces to that of calculating (exp jrx) where, from (86) and (96), 



= f c(x)dx. (102) 



We have seen in Section IV how to carry out this calculation if c(x) is 
a Gaussian or a random telegraph process. Similar remarks apply also 
to the calculation of the expected transmitted and reflected powers, 
and their variances. We do not give the results for the general case, 
although the calculations are straightforward, since the final expressions 
are somewhat lengthy. 

However, we will consider the case of matched lines, for which no 
reflections occur. Thus, with n P = 1, i> p = 1, (p = 1, 2), we have, from 
(97) through (99), 

T, =- e' Tl cos x, T 2 = JKe- Tl sin X - (103) 

This is for unit input voltage on line 1. Interchanging subscripts, for 
unit input voltage on line 2 we have 

T 1 = I e~ ri sin x , T 2 = e~ ri cos x , (104) 

K 

using (96). Thus, if V^O) = v x and ^(O) = v 2 , then 

Vl (l) = 7\ = e-* 1 ^ cos x + ^ 2 sin x ) , 

V 2 (l) = T 2 = e' ''(>•?-, sin x + v 2 cos x). (105) 



Note that 



T t | L " + | T 2 f = e - <r + 1 ""(| KVl | 2 + | v 2 \ 2 ). (106) 
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We let 

c(x) = cN(x), (107) 

where N(x) is a dimensionless stochastic process with zero mean. Then, 
from (37) and (102), 

x = c6(l). (108) 

Let us consider the case when N(x) is a Gaussian process. Then, from 
(52), 

(e" x ) = e-*'""' 1 ", (109) 

where <r 2 (l) is given by (51). Hence, from (105), 

{Tl ) = M - r V* , " ,c, \ (T 2 ) = v* vl e-*" Hl) . (110) 

Also, after some algebra, it is found that 

<l r, l 2 > = eJ ^Z [(I «, I 2 + I v* I 2 ) + (I *« I 2 - I * IV"" 1 "]. 
2 1*1 (111) 

(I r 2 I 2 ) = K ,r+r,, '[(l^ I 2 + I* I 2 ) - (1*1 I 2 - \ v > IV"""*]. 

(112) 
and 
<|«T, |*>-<|«T, I 2 ) 2 = <l r 2 I 4 ) - <l ^ IT 

= fcr ,<r+M, [i -«- 4,vc,) ]f|«C-«W | 2 [i + e - 4c "" ( "] 

-2(| y2 | 2 - |* | a )V 4o '" (,) }. (113) 

The first equality in (113) is a consequence of (106). From (111) through 
(113) it follows that the ratio of the standard deviation to the mean of 
\T P | 2 approaches 

1 2 22 | 

as c"V(Z) — > °° , both for p = 1 and ji = 2. 

Analogous results may be obtained when AT(x) is the random telegraph 
process, by using (73), and also in the case of weak general coupling and 
long sections of the coupled lines, by using (82). We remark that we have 
previously 2 calculated the average modal powers in these cases for two 
coupled lines with unequal propagation constants, corresponding to 
equation (16), using entirely different methods. 
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VI. TIME DOMAIN STATISTICS 



In this section we conclude the study of our models by calculating 
some of their time domain statistics. Consider first the single line and an 
incident wave moving to the right in x < of the form 

V(t - 7( ,r) = -7= r V(w)e ial, - y ° z) c/co, 

V 27T J-=o 

I(t - To-r) = ~ V(t - 7o.r). (115) 

Then the transmitted wave in x > I is 

v T (t - y,x) = -^ f^ T(co, z)i;(u>r (, ^' (I - n) C/co, 

Ir(t - 71 X) = jr V r (t - y lX ), (116) 

where T(co, /) is the transmission coefficient. If we substitute expression 
(41) for T(u, I) into (116) and formally interchange summation and 
integration, we obtain 

V T (t - y t x) = xf]p r "4jF= f V(o>) exp [ju[t - 7l (x - I) 

r = V 27T ■'-oo 

- (2r+ l)y(l + ed(l))]\ dco 



= X E pV(t - y,(.r - - (2r + 1) 7 (/ + e*(0)). 

(117) 



Therefore, 



<F T (< - 7 /.r)> = X £ p r (T(/ - yi (x - /) - (2r + l)y(l + «0(0))>. 

(118) 

In this formulation, the randomness appears just as we should expect in 
a random change in the electrical length of the central transmission line. 
From (117) it follows that 



(I V T (t - 7l x) |") 



= X" 



£ p r V(t - 7,(.r - - (2r + l)y(l + ed(l))) 



(119) 



To better understand some of the implications of these formulas, 
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let V(u) be a rectangular pulse, 

V(u) = 1, ^ U ^ T, 

= 0, u < 0, t < u. (120) 

Then if N(x) is one of the Gaussian processes discussed in Section IV, 
(V(t - 7l (x - I) - (2r + l)y(l + e0(Z)))> 

L f e - (B ' /2a,) V(t - yi(x -I)- (2r + l)y(l + e(f)) dd 

aV2ir J— 

i / ,I(-7 1 (i-l)-(2r + l)7'/^2(2r + l)«T»(01 

= J_ / _ e~ u ' dw 

V7T «'[«-7i(i-0-(2r + l)7l-'-/v'2(2r + l)«7»(J)l 



- \ [ erfc 



* - 7»(x - I) ~ (2r + W - t 



V2" (2r + 1)^(0 

_ erfc U - 7,(» ~ I) - (fr + Ml . (121) 

ertC l V§ (2r + l)ey<r(0 J J 

In (121) erfc (x) is the complementary error function." 

Equation (121) shows that the average field violates causality, since 
at x = I, for example, (V(t - (2r + l)y(l + ^(l)))) is positive for all 
— oo < t < oo . This is really a consequence of the fact that for a Gaus- 
sian process, at any point x, there is a positive probability that a sample 
function is less than -e" 1 . Hence, in our model the inductance and 
capacitance can both become negative, leading to the violation of 
causality. 

If N(x) is the random telegraph process, then 6(1) has the probability 
density function given in (71). It follows that in this case 

(V(t - 7l (x - I) - (2r + \)y(l + e0(Z)))> 

= \e- bl \v{t - 7l (x -I)- (2r + 1) 7 Z(1 + e)) 
+ V(t - 7<(z - I) ~ (2r + 1M1 - *)) 
+ b /' \l [b(? - 2 )*] 

+ Ui tf l i~e4 )i] } nt - y,{x ~ l) ~ (2r + 1)y(l + e6)) d6 ~\ ' 

(122) 
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If V(u) is the rectangular pulse (120) which arrives at x = at t = 0, 
then from (118) and (122), the average transmitted signal is a train of 
nonrectangular pulses. It is seen easily from (122) that 

(V(t - y t (x - I) - (2r + l)y(l + e0(O))> 

is a pulse which begins at t = y,(x - I) + (2r + 1)7/(1 - c) and ends at 
t = y,(x - I) + (2r + 1)7/(1 + e) + t. This shows that as long as 
e < 1, causality is preserved for the average transmitted signal since 
then 7,(x - I) + yl(l _ e ) > for all I ^ x. 

The duration of the rth average transmitted pulse is At = r -j- 
2(2r + \)yU, while the time between the end of the rth pulse and 
beginning of the (r + l)th pulse is 2?Z - 4(r + \)yU - T . Thus, each 
pulse in the train is longer than the pulse which preceeded it and all of 
these are longer than the incident pulse. Furthermore, no matter how 
short the incident pulse is, the average transmitted pulses eventually 
begin to overlap. 

The transmitted power can be treated in the same way. However, 
due to the overlapping of the pulses, the analysis is tedious and we do not 
discuss it here. 

Analogous results may be obtained in the case of two coupled trans- 
mission lines with identical propagation constants, if we consider 
lossless lines and suppose that the phase constants and the coupling 
coefficient are proportional to the frequency, so that 

T = juy, c = wye, (123) 

and 

r p0 = jcoy P o , T pl = juy pl , (p - 1, 2). (124) 

It is also assumed that the characteristic impedances K v , K p0 , and 
Kpi > (P = 1»2) are independent of co. We will confine our attention to 
the case of matched lines, so that no reflected waves occur. 

Let us consider an incident wave on line 1 moving to the right in 
x < 0, and of the form 

V l {t - 7io.r) ■■ ~ f° 'U 1 (a,y u( '- r ■ 0I, rfco, 

V 2tt J-* 

Hi - 7,o.r) = -£- V& - 7l0 .r). (125) 

-"-10 

Then, from (93), (94), (96), (103), (108), and (123) through (125), the 
transmitted waves in x > I are 



V2^ 
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V 1T (t - yux) = K u Iir(t - Yuz) 

= -±= P 1) 1 (co)e-'' U7Z co8(co 7 e^)y u( '- 7,,<1 "' >) &» 
V2tt -'-oo 

= h[Vt(t - y u (x - I) - uy[l - ed(l)]) 

+ Vt(t ~ yu(x - I) - uy[l + e6(l)]}, (126) 

and 

V 2T {t — y 2 ix) = K 2l I 2T (t — y2ix) 

h= (¥l f VxW 1mjl b^ («Y^(QV w(, " 7,,<8 " n) ^ 

'27T ^m/ ■'-«" 

= | (§j)Vi(' - ?«(* - - «r[l - tem 

- V 1 (t - y 2l (x -I)- uy[l + €0(0])} . (127) 

The transmitted waves in x > I corresponding to an incident wave 
on line 2 moving to the right in x < 0, and of the form 

V 2 (t - y 20 x) = -7= f^V 2 (co)e iai, - y " x) da, 

hit - y 20 x) = ^- V 2 (t - y 20 x), (128) 

are obtained by interchanging the subscripts 1 and 2 in (126) and (127). 
The transmitted waves corresponding to incident waves on both lines 
are obtained by linear superposition. The ensemble averages of the 
transmitted waves may be calculated as before. 

APPENDIX 

We here apply a limit theorem due to R. Z. Khas'minskii, 3 in order to 
determine the limiting probability density function for the process 
60(£/e 2 ), for bounded £, and e -» 0. Now, from (37), 

4- [tO{x)] = eN(x), e0(O) = 0. (129) 

ax 

We assume that N(x) is a bounded, zero-mean, wide-sense stationary 
stochastic process, with correlation function given by (58). Since is 
a scalar rather than a vector, and moreover the derivative of depends 
only on x, and not on 0, we have just about the simplest nontrivial 
application of the limit theorem. 
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There is only one quantity to be considered, namely, 

a(x, y) = (N{x)N{y)) = g(\ x - y |), (130) 

from (58). Then, in accordance with Khas'minskii's definition, 

a = lim \x J J a ^ X ' ^ dx dy 

= 2 lim [^ £ (X - z)g(z) &J , (131) 



X-oo 



after some integrations by parts. As required by the hypotheses of the 
limit theorem, a is independent of x . It is also required that the 
stochastic process N(x) satisfy a certain strong mixing condition, and 
the reader is referred to Khas'minskii's paper for a precise statement 
of this condition. 

For the case under consideration, the limit theorem states that, on 
the interval ^ £ ^ £ , where £„ is an arbitrary positive number, the 
process €0(£/e 2 ) converges weakly as e — » to a Markov diffusion process 
®(£) with zero drift and diffusion coefficient a. The drift coefficient R is 
zero since the right-hand side of eq. (129) is independent of 0. The 
probability density function of the limit process satisfies the equation 

% = \ a d<t 2 ' p(©>o) = a(e), (132) 

in view of the initial condition 0(0) = 0. Thus, 
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